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RESUMEN 

La posibilidad de que las erupciones de rayos gamma (GRBs) no sean emisiones isotropicas fue considerada 
teoricamcnte para atenuar el problema asociado a la gran cantidad de energfa implicada por el modelo de bola 
de fuego estandar para estos potentes fenomcnos. Sin embargo, el mecanismo por el cual, tras la deposion cuasi- 
isotropica de unos pocos 10 50 erg se origina una eyeccion colimada de plasma no pudo ser explicada de forma 
satisfactoria anah'ticamcnte. La razon de ello radica en que la colimacion dc un flujo salicnte por su sistema 
progenitor depende de una dinamica no lineal muy complcja. Ello ha hecho nccesario el uso de simulacioncs 
numcricas para arrojar algo de luz sobre la viabilidad dc algunos de los progenitores mas probables dc GRBs. 
En esta contribution rcvisare los hechos mas relevantes mostrados por tales simulacioncs numcricas y como 
estas han sido utilizadas para validar el modelos de estrella eolapsante (para GRBs largos) y el modelo que 
implica la fusion dc un sistema binario de objetos compactos (para GRBs cortos). 

ABSTRACT 

The possibility that gamma-ray bursts (GRBs) were not isotropic emissions was devised theoretically as a way 
to ameliorate the huge energetic budget implied by the standard fireball model for these powerful phenomena. 
However, the mechanism by which after the quasy-isotropic release of a few 10 50 erg yields a collimatcd ejection 
of plasma could not be satisfactory explained analytically. The reason being that the collimation of an outflow 
by its progenitor system depends on a very complex and non- linear dynamics. That has made necessary the 
use of numerical simulations in order to shed some light on the viability of some likely progenitors of GRBs. 
In this contribution I will review the most relevant features shown by these numerical simulations and how 
they have been used to validate the collapsar model (for long GRBs) and the model involving the merger of 
compact binaries (for short GRBs). 
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1. INTRODUCTION 
Our current understanding is that gamma-ray 
bursts (GRBs) are produced in the course of the 
birth of a stellar-mass black hole (BH). In other as- 
trophysical systems where accreting BHs fuel col- 
limated beams of plasma (e.g., AGNs and BH-X- 
ray binaries), there is a direct evidence for relativis- 
tic outflows and jet collimation which comes from 
the imaging of the system. Therefore, it seems 
reasonable to assume as starting point, that also 
GRBs result from relativistic, collimated outflows 
from accreting, stellar-mass BHs. We have inferred 
that outflows yielding GRBs are relativistic because 
of a couple of observational constraints, namely, 
the detection of radio scintillation of the intcrstel- 



cxpansion is also necessary to overcome the theoret- 
i cal constrain imposed by the compactness problem 
(jCavallo fc Reesl . Il978j ). However, we only have an 



indirect evidence of collimation based on the obser- 
vational constraint posed by the achromatic break 
in the afterglow ligh t curve of some GRBs (e.g., 



Harrison et al. . 19991 ). From the theoretical point of 
view, if GRBs are collimated events, the true emit- 
ted energy is reduced by a factor /n — 2 /2 



(|RhoadsLll999l : ISari et al.l . ll999f ). i.e., £ 7 = faE. 



lar medium ( Frail et all 119970 and the measure 



ment of s uperluminal proper motions in imaged af- 



terglows (j Tavlor et all |200J). The ultrarclativistic 
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where E~ r - lso is the detected equivalent isotropic en- 
ergy. Nonetheless, the mechanism by which after 
the quasi-isotropic release of an amount of energy 
in the range 10 48 — 10 51 erg results in a collimated 
ejection has not satisfactory been explained analyt- 
ically. The reason being that the collimation of an 
outflow by the progenitor system depends on a very 
complex and non-linear dynamics. That has made 
necessary the use of numerical simulations in order 
to understand the collimation mechanism as well as 
to shed some light on the viability of some systems 
proposed to be the progenitors of GRBs. A robust 
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result obtained in numerical simulations of genera- 
tion of GRBs is that the progenitor system yields 
collimated outflows under rather general conditions 
independen t on whether the outflow is initiated ther- 
mally (e.g., Alov et al. . 2000l 2005f ) or it is magnet- 
ically driven ( McKinnevT 2006 ) . 

In this contribution I will review the most rele- 
vant features shown by these numerical simulations 
and how they have been used to validate the col- 
lapsar model (for long GRBs; § [2]) and the model 
involving the merger of compact binaries (for short 
GRBs; §E]). 

2. OUTFLOWS EMERGING FROM 
PROGENITORS OF LONG GRBS 

Among the plethora of models devised to ex- 
plain the origin of long GRBs (lGRBs), the most 
widely accepted was put forward by Woosley. In the 
original coll apsar model, als o known as failed super- 
nova model (|wooslevl . [l993h . the collapse of a mas- 
sive (Mzams ~ 30 M Q ) rotating star that does not 
form a successful supernova but collapses to a BH 
(Mbh ~ 3 M©) surrounded by a thick disk. The vis- 
cous accretion of the disk matter onto the BH yields 
a strong heating that, in its turn, produces a co- 
pious amount of thermal neutrinos and antineutri- 
nos, which annihilate preferentially around the ro- 
tation axis producing a fireball of e + e~ pairs and 
high energy photons. Later it was noted that, per- 
haps z/-powered fireballs might not be sufficiently 
energetic to fuel the most powerful GRB events 
and, thus, the collapsar model was extended to ac- 
count for alternative ene rgy extraction mechanisms 



(jMacFadven et all l200lh . More cxplicitcly, the ac- 



cretion energy of the torus could be tap by suffi- 
ciently strong magnetic fields (hydromagnetic gen- 
eration) by means of the Blandford-Payne process 
dBlandford fc Pavnd . ll982l ). or a non-trivial fraction 
of the rotational energy of the BH may also be con- 
verte d into a Poynting flux (jBlandford fc Znajekl . 
1977h . 

The scape of the newly born fireball and its ter- 
minal Lorentz factor (Too) depend on structural and 
on dynamical factors. The critical structural fac- 
tors are the environmental baryon density in the 
funnel around the rotation axis of the star and 
the ability of the progenitor star to loose its outer 
Hydrogen envelope. An under-dense funnel forms 
along with the accretion torus if the specific an- 
gular momentum of the core of the star lies in 
the range 3 x lO^cm 2 ^" 1 £ j < 2 x 10 17 cm 2 s" 1 
( MacFadven fc Wooslevl . [l999h . The existence of the 
funnel is key to collimate the fireball and to permit 



its propagation through the progenitor. The most 
favourable conditions for the generation and prop- 
agation of the fireball happen when the density of 
the funnel (pf) is much smaller than the density of 
the torus (ptorus), namely, Pf/ ptoms ~ 10~ 4 - 10~ 3 
dMacFadven fc Wooslevl . 1 1993) . The likelihood that 



the progenitor star had lost its Hydrogen envelope 
depends on a number of factors which occurrence 
is still a matter of debate like, e.g., the genera- 
tion of stellar winds, the interact ion with a com- 
panion (jPodsiadlowski et al I l2004i etc. The rele- 
vance of the lost of the hydrogen envelope resides 
on the fact that, unless the densi ty of the funnel is 
extre mely small (as proposed by IMeszaros fc Reed . 



20011 ). only a mildly relativistic, poorly-collimatcd 



fireball would reach the outer edge of the hydro- 
gen envelop after very long times and with relatively 
small Lorentz factors (r ~ 2) implying that the ob- 
servational signature would be an X-ray/UV tran- 
sient with a duration of ~ 10 — 1000 s but not a 
GRB (jMacFadven et all . 1200 ll ). Finally, the most 



important dynamical factor setting Too is the the 
amount of baryons entrained as the fireball propa- 
gates through the stellar mantle. 

Several analytic works have made estimates 
about the collimation angle and the Lorentz fac- 
tor when the fireball breaks out the surface of the 
star (e.g. 



IMeszaros fc Reesl . 120011 ). Nevertheless, 
the complexity inherent to the non-linear (mag- 
neto)hydrodynamic interaction of the fireball plasma 
with the stellar environment makes it unavoidable 
the use of numerical simulations. With these simu- 
lations we have been able to give preliminary answers 
to the following questions: 

Collimation. The generated outflows are inertially 
or magnetically confined. In the first case, the colli- 
mator is the funnel within the progenitor while in 
the second case, the flow is sclf-collimatcd by its 
own magnetic field if it is strong enough. Rather 
independently on the initial conditions and on the 
inclusion of magnetic fields, the typical outflow half- 
opening angles, when the jet reaches the surface of 
the progenitor star, are $brcak ~ 5°. These small 
half-opening angles result from the recollimation of 
the outflow within the progenitor and they are in- 
dependent on whether the boundary conditions are 
set to initiate the jet wit h much large r half-o pening 
angles (e.g., 6 = 20°; IZhang et all 120031 ) or on 
whether the jet is generated by an energy release 
into a volume spanning a half-opening angle 9d = 30° 



much larger than #t>reak ( Alov et al. . 2000l ). In the 



course of their propagation through the progenitor 
the jets develop a non-homogeneous structure, trans- 
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verse to the direction of motion, whose main fea- 
tures are an internal ultrarelativistic spine (where 
the Lorentz factor may reach r coro ~ 30 — 50 at jet 
breakout) within a half-opening angle of < 2° later- 
ally endowed by a moderately relativistic, hot shear 
layer (T sh i ~ 5-1 0) extending up to 6» shl < 20° - 30° 
(|Alov et al.l . l2002f ). We point out that the transverse 
structure of the jet is nearly Gaussian both in sim- 
ulations including magnetic field s (e.g., McKinne y 



2006) or not including them (e.g jAlov et alll2000f) 



However, a more accurate fit of the transverse struc- 
ture of the jet, c annot be ac commodated by a simple 
Gauss function jAloyj, 1200 ll ) . 

Variability. All produced outflows are highly 
variable due to the generation of Ke l vin-H elmholtz 
(lAlov et all. l2000l: [ Gomez fc Harded . 12004 ) , shear- 
driven (|Alov et all . 20021 ) or pinch magnetohy drody- 
namic (MHD) instabilities (jMcKinnevi . 12006( 1. Such 
extrinsic variability is independent on the [intrin- 
sic) variability of the energy source and leads to 
the formation of irregularities in the flow which are 
the seeds of internal shocks in the outflow. Ex- 
cept in cases in which the source may produce 
quasi-periodic variability (perhaps induced by pre- 
cession or nutation modes of the accretion disc), 
the extrinsic variability might be indistinguishable 
form the intrinsic one. Numerical simulations of 
three-dimensional (3D) relativistic jets propagating 
through collapsar-l ike environments s how that such 
jets are also stable (jZhang et all 120041 ) but it still re- 
mains unknown whether 3D relativistic, magnctohy- 
drodynamic, collapsar-jets will also be stable along 
its whole trajectory. 

Jet breakout. The jets generated are much lighter 
than their baryon reach environments. Thus, they 
propagate through the collapsar at moderate speeds 
(~ 0.3c) and fill up thick cigar-shaped cavities or co- 
coons of shocked matter that also propagates along 
with the beam of the jet and that, eventually may 
break the surface of the collapsar. Since in the 
cocoon a few 10 50 erg may be stored as the jet 
drills its way through the star, it has been pro- 
posed that its eruption through the collapsar surface 
could yield a number of 7-r ay/X-ray/UV-transients 



( Ramirez- Ruiz et al. . 20021 ). Indeed, it has been 



proposed that GRBs are but one observable phe- 
nomenon accompanying black hole birth and other 
possibilities may arise depending on the observer's 
viewing angle with res pect to the pr opagation of 
the ultrarelativistic jet ( Wooslev . 2000h . Thus, in a 
sort of unification theory for high-energy transients, 
one may see progressively softer events ranging from 
GRBs (when the jet emergence is seen almost head 



on) to UV flashes (when the jet eruption is seen 
at relatively large polar angles) and accounting for 
X-ray rich GRBs (XRR-GRBs) and X-ray flashes 
(XRF) at intermediate angles. The jet emergence 
through the stellar surface and its interaction with 
the stellar wind (which likely happens during the 
late stages of the evolution of m assive stars) could 



lead to some precursor activity (jMacFadven et al 
20011 ). Furthermore, ^-powered jets are very hot at 
breakout (~ 80% of the total energy is stored in the 
form of thermal energy) which implies that jets can 
still experience an additional acceleration by conver- 
sion of thermal into kinetic energy, even if the energy 
source has ceased its activity. 

2.1. IGRBs produced in collapsars: MHD- or 
neutrino-powered jets ? 

From a aesthetic point of view, it is beautiful 
if any invoked jet powering mechanism explains not 
only events with relatively small Lorentz factors T ~ 
100 but also the occurrence of events with very large 
inferred values of T ~ 500 (|Lithwick fc Saril . l200lh 
or eve n T ~ 1000 as suggested by mode ls for some 
GRBs (jSoderberg fc Ramirez-Ruizl . I2003T ) . However, 
there are no fundamental reasons to argue in favour 
of a unique and universal mechanism to extract the 
energy stored in the progenitor. 

Purely hydrodynamic, i^-powered jets in collap- 
sars of type-I seem to be able only to produce more 
moderate terminal v alues of the bulk Lorentz facto r 
(r ~ 100 - 400; e.g., lAlov et all l200d : IZhang et al 



20031 ). even if there is a further acceleration of the 
forward shock as a result of an appropriate den- 
sity gra dient in the mediu m surrounding the pro- 
genitor ( Alov et al. . 2000h . unless the density in 
the funnel around the r otation axis is very small 
( Meszaros fc Reed . l200lh . In collapsars of typc-II 
the mass accretion rate [M ~ 10 -3 — 10 -2 M Q s _1 ) 
is insufficient to produce a i^-powered jet, but it 
may suffice to generate MHD-powered jets whose 
observational signature might be weak, poorly col- 
limated X-ray/UV transients or very long (~ 100 — 
1000 s) GRBs of low E iso and small L if the pro- 
genitor star is able to lo ose its Hydrogen envelope 



(jMacFadven et all 120011 ). 



From extrapolations of the numerical results 
of axisymmetric j ets generated electromagnctically 
(McKinneyJ, 120061 ) very large values of the asymp- 



totic Lorentz factor of the outflow (r^ ~ 1000) can 
be attained. Thus, it seems rather plausible that 
MHD mechanisms have to be employed to generate 
jets with T - 500 - 1000 (jMcKinnevi |2005| . l2006h 
and, invoking the existence of a universal energy ex- 
traction mechanism, one may also argue that also 
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Fig. 1. Logarithm of the density (gray-scale) with overimposed magnetic field lines (white lines) and total magnetic 
field strength (black contours) corresponding to models collapsing stellar cores with a small initial rotational energy and 
rotating almost rigidly (left) and with a larger initial rotational energy and differentially rotating (right). Each panel 
displays the state of the stellar core after the bounce for four different initial magnetic field strengths (Bo = 10 10 G, 
10 11 G, fO 12 G and 10 13 G in the clockwise direction starting from the upper left corner). The figures correspond to 
models of Obergaulinger et al. (2006a). 



jets with much more moderate terminal Lorentz fac- 
tors are also MHD-generated. Nevertheless, these 
arguments have a number of issues: 

1. - The estimates of the terminal Lorentz factor 
of MHD-generated jets are based on axisymmctric 
models. The 3D-stability of relativistic magnetized 
jets is still a matter of debate which needs of much 
more complex numerical simulations than the ones 
produced so far. Should 3D-MHD jets be unsta- 
ble, should the terminal Lorentz factor not reliably 
be calculated with the state-of-the-art axisymmetric, 
general relativistic, MHD (GRMHD) simulations up 
to date (|McKinnevl . l2006t) . 

2. - GRMHD initial models of accreting BH sys- 
tems consist on more or less realistic matter distri- 
butions over which an assumed poloidal field is im- 
posed, i.e., they are not the result of a consistent 



magne to - rotational core c ollapse (e.g.. lMizuno et al 



2004al lbl: iMcKinneyi I2006T ) . The use of poloidal fields 
is triggered by the fact that purely toroidal field con- 
figur ations do not yie l d the production of bipolar 
jets (|De Villiers et all . l2005f) . On the other hand, 



the initial magnetic field strengths are assumed, not 
consistently computed from the core collapse of mas- 
sive stars. Typically, initial field strengths as large 
as B ~ 10 15 — 10 16 G are used. Such large val- 
ues of B in combination with maximal values of 



the dimensionless angular momentum of the BH 
C?bh ~ 1) necessary in order to efficiently ex- 
tract energy via Blandford-Znajek (BZ) mechanism, 
because the BZ-power scales rather sensitive with B 
as £ R7 ~ 10 50 7g „(M R H/3M^) 2 (g/10 15 G) 2 - 



(Lee et al.. 200 



/RH 



' ergs 



3.- When numerical simulations of the magne- 
torotational core collapse of massive stars ar e per- 
formed (see, e.g.. lObergaulinger et al. . 2006bl la. and 
references therein), bipolar outflows as well as thick 
toroidal structures surrounding a central, low density 
funnel in the collapsed core are only generated for 
initial magnetic field strengths Bo ^ 10 12 G and for 
rather fasl|f| and differentially rotating stellar cores 
(Fig.[T]right). Such initial field strengths and angular 
velocities are well beyond the ones predicted by the 
state-of-the-art calc ulations of rotating massive stars 
dHeger et all 120051) . These initially strongly magne- 
tized models likely develop collapsed cores, with a 
mass of M c ~ 1 Mq and a specific angular momen- 
tum J c r~j 10 cm s , and may form a rapidly ro- 
tating BH with j'bh ~ 1. Furthermore, the winding 
up of the initial poloidal field leads to maximal field 
strengths B max ~ 10 15 G, being the field predomi- 
nantly toroidal Storoidai/-B P oioidai ~ 1 — 10. On the 



3 The rotational energy being as large as ~ 4% of the grav- 
itational energy. 
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other hand, initially rigid, more moderately rotating 
cores do not yield tori around low-density funnels, 
do not produce bipolar jets and the maximum field 
strengths are ~5x 10 14 G (Fig. ED left). Considering 
that the collapsed cores are smaller (M c ~ 0.75 Mq) 
and slowly rotating (J c ~2x 10 15 cm 2 s _1 ) than in 
the previous case one may expect that the newly 
born BH resulting from the posterior evolution of 
these kind of cores will not be a maximally rotating 
but, instead, they will form BHs with more moder- 
ate j'bh ~ 0.6. Thereby, the conclusion that seems to 
emerge from detailed numerical simulations of Ober- 
gaulinger et al. is that if the initial magnetic field 
strength and rotational energy is as small as pre- 
dicted by the most detailed stellar evolution models, 
the collapsed core does not hold the appropriate con- 
ditions (B 10 15 G, jbh ~ 1) to efficiently extract 
energy via BZ-mechanism and, conversely, unrealis- 
tically large initial magnetic fields and rotational en- 
ergies need to be invoked to expect a BZ-like mech- 
anism to operate efficiently. However, we have to 
be cautious in order not to extract too far reaching 
conclusions from the previously mentioned numer- 
ical work. It remains true that the results of the 
simulations of Obergaulinger et al. are handicapped 
because they do not include general relativity, be- 
cause the numerical resolution is still small to cap- 
ture all the relevant magneto-rotationally unstable 
modes and, because they are restricted to axisym- 
metric models. 

From the above points, one may infer that the 
dynamical relevance of the magnetic field in the pro- 
cess of energy extraction from the central source 
will depend on fine details of the magnetorotational 
collapse of the collapsar core. On the other hand, 
the process of vP-annihilation as the primary source 
of energy that fuels an ultrarelativistic fireball also 
needs of a more careful study in order to know how 
much energy such a process may release in the pro- 
genitor system and how such an amount of energy 
depends on the physical conditions of the progenitor. 

A step towards such goal is the work of Birkl et 
al. (2006), which contributes to better understand 
how the energy deposition rate due to the process of 
^^-annihilation (E v p) depends on general rclativis- 
tic (GR) effects and on different neutrinosphere ge- 
ometries in hyperaccreting stellar-mass BH systems. 
Birkl et al. consider two families of neutrinosphcrcs. 
On the one side, idealized geometries as thin disks, 
tori, and spheres. On the other side, more realistic 
models are constructed as non-selfgravitating equi- 
librium matter distributions for varied BH rotation. 
Independent of whether GR effects are included, con- 



sidering the same values of temperature and sur- 
face area for an isothermal neutrinosphere, thin disk 
models yield the highest energy deposition rates by 
^-annihilation, while spherical neutrinospheres lead 
to the lowest ones. Considering isothermal neutri- 
nospheres with the same temperature and surface 
area, it turns out that compared to Newtonian cal- 
culations, GR effects increase the annihilation rate 
measured by an observer at infinity by a factor of 
2 when the neutrinosph ere is a disk (in agreemen t 
with the previous works; lAsano fc Fukuvamal . 120011 ). 
However, in case of a torus and a sphere the influ- 
ence of GR effects is globally only ~ 25%, although 
locally, particularly in the vicinity of the rotation 
axis of the system, it can be significantly larger. Fo- 
cusing on the dependence of the energy deposition 
rate on the value of Jbh, it is found that increas- 
ing it from to 1 enhances the energy deposition 
rate measured by an observer at infinity by roughly 
a factor of 2 due to the change of the inner radius of 
the neutrinosphere. Furthermore, although the ab- 
solute values of the energy deposition rate have to be 
taken with care (because of the steady state approx- 
imation used and the need o f more real i stic m od- 
els for the accretion disk; see Birkl et al.l . I2006L for 
accretion disks of mass similar to the one expected 
in the collapsar model (Mdi s k ~ 0.01 Mq) typically, 
Evv ~ 10 50 — 10 51 ergs _1 . Even if only an small 
fraction (~ 10%) of that energy were used to boost 
a polar outflow, there is fair chance for neutrinos 
to be the dominant energy source of the fireball (at 
least, in some cases when the magnetic field is not 
too large). The most likely scenario that can be de- 
vised is that both mechanisms (MHD and neutrino 
energy release) might be operating simultaneously. 
Indeed, for MHD-produced jets, neutrinos will play 
a fundamental role in the p air-loading of the jet (e.g., 
Levinson fc Eichlerl . Il993f ) while, for ^-powered jets, 



the magnetic field may be important to collimatc the 
thermally generated outflow. Which of the two en- 
ergy deposition mechanisms dominates in every sin- 
gle GRB will depend on the exact conditions in the 
precollapse progenitor. 

3. OUTFLOWS EMERGING FROM 
PROGENITORS OF SHORT GRBS 

Nowadays, it is commonly believed that short 
GRBs (sGRBs) are generated after the merger of 
a system compact binaries formed by either two 
neutron stars (NS s) or a ne u tron s t ar and a BH 
( Paczvnskil. Il986t iGoodmanl. Il98& lEichler et af 



19891 ; iMochkovitch et all . Il993l ). The remnant left 



by the merger consists of a newly born BH black 
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hole girded by a thick gas torus from which it swal- 
lows matter at a hypercritical rate. In such situa- 
tion radiation is advected inward with the accretion 
flow and the cooling is dom i nated by the emission 



of neutrinos ( Popham et al. . 19991 ). As in the case 



of progenitors of lGRBs, these neutrinos might ei- 
ther be the primary energy source blowing a fire- 
ball of e + e~ pairs and photons or to act as media- 
tor in hydromagnetic or electromagnetic energy ex- 
traction mechanisms (see Sect. [5]) to pair-load the 
Poynting dominated outflow. A fundamental dif- 
ference with respect to the collapsar model is that 
the accretion thick disk is cannot be continuously 
refilled from a surrounding matter reservoir (the stel- 
lar matter in case of a collapsar) and, therefore, the 
duration o f the produced outflows is, in part (sec 
Alov et all fiooBh roughly limited by the time dur- 
ing needed by the black hole to engulf most of the 
matter of the accretion disk, namely, a few 100 ms. 
This limit on the time scale, set by the ON time t ce 
of the source, of holds for both ^-powered jets and 
for MHD-generated outflows. In the first case, the 
neutrino luminosity fades as the mass of the disk de- 
creases and, thereby, there will be a critical torus 
mass below which a plasma outflow cannot be sus- 
tained. In the second case, for the same reason, 
there will not be sufficient neutrinos that pair-load 
the Poynting dominated jet after a sizable fraction 
of the torus has been accreted. 



10 erg above the poles of a stellar mass BH 
in a region of nearly vacuum may yield an ul- 
trarelativistic outflow, numerical simulations are 
needed to attempt to answer questions about the 
collimation mechanism of the polar outflow, the 
opening angle of the ultrarelativistic ejecta, the 
asymptotic Lorentz factors that can be attained, 
the internal structure of the outflow, the dura- 
tion of a possible GRB event, and the isotropic 
equivalent energy which an observer would infer 
by assuming the source to expand isotropically. 
These questions may at most be guessed, but they 
cannot be reliably answered on grounds of merger 
models and a consideration of their energy release 
by neutrino emission and the subsequent conversion 
of som e of this energy by ^- a nnihilation t o e + e- 
pairs dRuffert fc Jankal. f 19991 Ijanka et all Il99£ : 



Rosswog fc Ramirez- Ruiz 



2003: Birkletal 



.20021 iRosswog et al 
20061 ). Only self-consistently 
time-dependent (magncto)hydrodynamic modeling 
may give us some insight on the former questions. 
The reason being that the relativistic outflow 
develops in a complex interaction with the accretion 



torus, cleaning its own axial funnel such that 
later energy deposition encounters a much reduced 
baryon pollution. 

Some keys to answer the questions mentioned in 
the previous paragraph have been very recently re- 
vealed by time-dependent numerical simulations in 
which the main results are: 

Collimation. The generated outflows which may 
yield GRB sig natures are either collimated by the 
accretion disk ( Alov et al. . 20051) or self-collimated 



by the magnetic field (jMcKinnevl . 120061) . depending 
on whether the jet is initiated thermally or magnet- 
ically, respectively. The typical outflow half-opening 
angles are ~ 3° — 25° and, as in the case of jets pro- 
duced in collapsars, the baryon-poor outflows dis- 
play a transverse structure. This structure shows 
a central core which spans a half opening angle 
0core < 3° - 12° where r corc > 100 flanked later- 
ally by a layer, extending up to # s hi ~ 25° where the 
Lorentz factor smoothly decays to moderately rel- 
ativistic values and where a sizable fraction of the 
total energy is stored. This layer is rather hot in 
thermally initiated outflows and has the potential of 
accelerating even after the energy release by the cen- 
tral engine has ceased. Similar to the jets produced 
in collapsar environments the transverse structure 
of the Lorentz factor could be roughly fit by Gaus- 
sian profiles. However, somewhat more complicated 
functions are requ ired to provide more accurate fits 



Although it seems likely that releasing a few (jAloy et al. |, |2005J) 



Variability. Even injecting energy close to the BH 
even horizon at constant rates, the produced out- 
flows are highly variable. The interaction of the 
newborn fireball with the acc retion torus yields the 
growth of Kelvin-Helmholtz ( Alov et al. . 20051 ) in- 



stabilities. The variability in case of MHD jets i s im- 
printed by pinch instabilities ( McKinnev . 20061 ). Up 
to date only axisymmetric models have been com- 
puted. All of which seem to be stable or marginally 
stable. It is not yet numerically verified whether 3D 
jets emerging from hyperaccreting BHs are stable. 
Influence of the environment. Mergers of com- 
pact objects may take usually place in the intergalac- 
tic medium or in the outer skirts of their host galax- 
ies. Thus, the environment of the merger may have 
a very low density. However, in the course of the 
merger, after the first contact of the two compact 
objects, there is an ejection of matter, which is larger 
close to the orbital plane. Such matter forms a cool 
baryon re ach environment with a mass Afhain < a few 
lfr_f M Q (|Ruffert fc Jankal . l200lHOechslin fc Jankal . 
20061 ). The exact amount of mass ejected sensitively 



depends on whether the two compact objects are 



RELATIVISTIC OUTFLOWS IN GRBS 



7 



both NSs or one of them is a BH, on the initial 
mass ratio between the two, etc. Thereby, when 
the central BH forms, the newly born system may, 
in some cases, be embedded into a relatively high 
density halo. Mergers in low density environments 
may fuel ultrarclativistic outflows with the poten- 
tial to produce normal sGRBs, while in case that 
the merger occurs in high density media, the obser- 
vational signature is not a sGRB but, most likely , 



a flas h in the soft X-ray or UV bands (jAlov et al 
2005). In the latter case, the resulting event might 
only be observable if it happens very close to the 
Milky Way. The fact that depending on the envi- 
ronmental density an sGRB can be produced or not 
has the direct implication that not every merger may 
yield a sGRB. This fact has to be considered when 
making estimates about the true rates of sGRBs and 
compared with the ra tes of NS+NS mergers (e.g., 
Guetta fc Pirar] . l2005h . 

Asymptotic Lorentz factor. The saturation value 
of the outflow Lorentz factor, r^, is difficult to es- 
timate on the basis of numerical simulations that 
cover only the initial fraction of a second in the evo- 
lution of an ultrarclativistic outflow generated in a 
hyperaccreting BH. Despite this difficulties rough es- 
timates can be made. For instan ce, in case that th e 
outflows are thermally generated ( Alov et al"1 . l2005l ). 



there is a clear trend to produce much higher val- 
ues of Too for sGRBs (r m > 500 - 1000) than for 
lGRBs (Fqo ~ 100). The reason for the difference re- 
sides on the much smaller density of the environment 
of the merger (even accounting for the mass ejected 
from the compact objects after the first contact; see 
above) as compared with the baryon-polluted envi- 
ronment that a relativistic jet finds inside a collaps- 
ing massive star. It has been speculated that this 
difference in Lorentz factor might be the reas on for 
the paucity of soft sGRBs (jjanka et all l2006h . The 
former trend seems not to be follow ed by MHD- 
generated outflows (|McKinnevl . l2006h . The most 
likely reason being that McKinney's simulations arc 
set to be scale free, while there should be a big dif- 
ference between the environments of mergers of com- 
pact objects and the interior of collapsars. Proba- 
bly, such a difference cannot be accommodated easily 
with simple scale- free power-laws for the distribution 
of the physical variables. 

Duration of the events. As pointed out by Aloy 
et al. (2005) and Janka et al. (2006), the "shells" 
ejected by the central engine, accelerate much faster 
in the leading part of the outflow than the shells in its 
lagging part. The rear shells therefore need a longer 
time to reach velocities v ~ c. This differential ac- 



celeration at early and late times of the relativistic 
jet leads to a stretching of the overall radial length 
of the outflow, A, relative to t cc times the speed of 
light c, A > ct cc . This stretching has the important 
consequence that the overall observable duration of 
the GRB (in the source frame), T = t& = A/c, may 
be a factor of 10 or more longer than t ce , even when 
the GRB is produced by internal shocks. 

4. SUMMARY 

The numerical modelling of progenitors of GRBs 
has allowed us to gain some insight into a num- 
ber of important issues related with the nature of 
the outflows produced by these systems. First, it 
has allowed us to verify that some (but probably 
not all) collapsars can yield collimated relativistic 
outflows that turn into lGRBs at large distances 
from the source. Likewise, only a fraction of the 
mergers of compact objects may yield sGRBs. Sec- 
ond, the numerical modeling has given us informa- 
tion about the collimation mechanism and typical 
outflow opening angles. Third, it has shown that 
the outflow is heterogeneous both along the direc- 
tion of propagation and transverse to it, even if the 
central engine releases energy at a constant rate. 
Fourth, it seems rather plausible that some lGRBs 
with very high Lorentz factor (r >> 100) need of a 
MHD jet-formation mechanism. On the other hand, 
some other lGRBs with more moderate Lorentz fac- 
tors (r ~ 100) can be explained by ^-powered jets, 
particularly if t cc < 10 s. Fifth, MHD- and v- 
mechanisms may work simultaneously. Therefore, 
it is likely that, in some cases MHD processes domi- 
nate the jet generation while in others neutrinos may 
be the dominant energy source. 

The numerical modeling done so far is still insu- 
ficient. In order to start from consistent initial mod- 
els, the state-of-the-art numerical codes must incor- 
porate, at least, the effects of strong gravity, mag- 
netic fields and a detailed neutrino transport. In the 
near future we may see how all these elements are in- 
cluded in more realistic numerical experiments that 
will deepen our understanding on how ultrarelativis- 
tic outflows are produced in progenitors of GRBs. 
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